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The Oide effect considers the synchrotron radiation in the final focnsing qnadrnpole and it sets a 
lower limit on the vertical beam size at the Interaction Point, particularly relevant for high energy 
linear colliders. The theory of the Oide effect was derived considering only the radiation in the 
focusing plane of the magnet. 

This article addresses the theoretical calculation of the radiation effect on the beam size consider¬ 
ing both focusing and defocusing planes of the quadrupole, refered to as 2D-Oide. The CLIC 3 TeV 
final quadrupole (QDO) and beam parameters are used to compare the theoretical results from the 
Oide effect and the 2D-Oide effect with particle tracking in PLACET. The 2D-oide demonstrates 
to be important as it increases by 17% the contribution to the beam size. 

Further insight into the aberrations induced by the synchrotron radiation opens the possibility to 
partially correct the 2D-Oide effect with octupole magnets. A beam size reduction of 4% is achieved 
in the simplest configuration, using a single octupole. 

PACS numbers: 29.20.Ej, 41.85.Lc 


I. INTRODUCTION 


with 


Synchrotron radiation in a focusing quadrupole mag¬ 
net of length L and gradient k, schematically represented 
in Fig. 1, changes the energy of the particle and modifies 
the focusing effect. This results in a limit on the mini¬ 
mum beam size at the Interaction Point (IP) located at 
a distance I* from the quadrupole. This is referred to as 
Oide effect [1]. 



FIG. 1. The solid and dashed lines represent the particle 
trajectory without and with radiation, respectively. 
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and Ae is the Compton wavelength of the electron, re 
is the classical electron radius and 7 is the relativistic 
factor. 

The primitive of the double integral used to calculate 
F in Eq. (2) is derived analytically in [2] with the goal 
to increase the computational calculation speed. It has 
been included in MAPCLASS2 [3-6] to be used in lattice 
design and optimization. 

Although the total contribution to beam size de¬ 
pends on the lattice and beam parameters, the minimum 
achievable beam size is given by [1] 


The beam size growth due to radiation is added 
quadratically to the linear beam size CTq = eP* where 
/3* represents the optical beta function at the IP and e 
is the emittance. Therefore, = a q + cr^j^g, where the 
beam size contribution from the Oide effect is [1], 
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(3) 

where cat = ye is the normalized emittance, showing the 
independence from beam energy. 

Table I shows relevant parameters of the last verti¬ 
cally focusing magnet QDO (fc, L), the beam and lattice 
optics (ew, 7, CTo, I*) for the two main linear collider 
projects ILC 500 GeV [7], and CLIC at 500 GeV and 
3 TeV [8]. These are used to show that the contribution 
of the Oide effect to vertical beam size is only significant 
for CLIC 3 TeV. Also, columns a and CTmin, show that 
the final vertical beam size is comparable to the minimum 
achievable. 

If none of the beam parameters or I* is to be changed 
then F can be used as a figure of merit of the quadrupole 
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CLIC 3 TeV 

20 

2935.0 

0.70 

0.116 

2.73 

3.5 

4.086 

0.87 

1.10 

1.00 

CLIC 500 GeV 

25 

489.2 

2.3 

0.077 

3.35 

4.3 

4.115 

0.08 

2.3 

1.17 

ILC 500 GeV 

40 

489.2 

5.7 

0.170 

2.20 

4.3 

9.567 

0.04 

5.7 

1.85 


TABLE I. Vertical beam size and radiation beam size contribution for three lattices, ejv is the normalized emittance, Civ = 7 £. 


set-up as it is calculated only from fc, L and I*, where the 
target is to reduce F as much as possible. The standard 
procedure of reducing F is by increasing the length of the 
quadrupole and reducing its gradient. In the following 
this is illustrated for CLIC 3 TeV. 

The relative increase in the beam size due to the Oide 
effect is given by 

= yi + f^oide/o-Q- (4) 

Figure 2(a) shows the the relative increase due to Oide 
effect when k is set to cancel the optics function alpha, 
a, just at the quadrupole opposite face to the IP, the 
particle input. This is done to set an absolute minimum 
in the Oide effect out of reach for realistic Final Doublet 
(FD) designs. Figure 2(b) shows the k values previously 
mentioned. 

The current QDO almost doubles the vertical beam 
size. It might be possible to reduce by a small fac¬ 
tor by increasing the current quadrupole length and pos¬ 
sibly using a lower k, however, this solution also points 
to increasing the lattice length. This could also generate 
difficulties in lattice design because of chromaticity and 
magnet stabilization. Quad lengths larger than 10 m do 
not lead to further significant improvements. 


II. 2D-OIDE EFFECT 

In order to gain understanding for the relevant case of 
CLIC at 3 TeV, Section II is dedicated to the theoretical 
derivation of the effect of radiation including the horizon¬ 
tal beam size and optics lattice parameters (2D-Oide), 
and comparing results with particle tracking. This leads 
to a possible way to mitigate the effect, alternative to 
enlarging QDO, consisting in removing the correlation at 
the IP via normal octupole magnets. This is tested for 
the simplest case of one octupole magnet as an example. 


A. Analytical derivation 

Being Eq the nominal energy of the beam and u the 
energy of the photon radiated, in a first order approxi¬ 
mation where u Eq, the total effect of radiation in the 
vertical displacement at the IP is calculated integrating 



(a) 



(b) 

FIG. 2. Oide effect beam size contribution for CLIC 3 TeV 
design parameters, (a) Total beam size normalized to de¬ 
signed linear beam size as a function of quad length for the 
minimum focusing k and the current QDO. (b) k in the two 
previous cases. 

along the magnet length as 

^y= [ [ fy{^si)^yQdsids, (5) 

Jo Jo ^0 

where the inner integral represents the kick given by the 
radiated photon at a certain location s propagated to 
the IP assuming a strong focusing in the vertical plane 
{I* » Py) [1]. y'o is the vertical particle angle at the IP 
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FIG. 3. The solid and dashed lines represent the particle 
trajectory without and with radiation, respectively. 


without radiation and 


fy {(j)) = sin 0 + -s/fcZ* cos (j). (6) 


Because u is much smaller than Eq in Eq. (5), it is still 
possible to use the result in [9] for the average energy loss 
per unit length in the magnet u and its second moment 
shown in Eqs. (7) and (8): 


u 2 7^ 
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This allows to calculate two important values: the mean 
effect of the radiation in the trajectory of one particle 
along the magnet Ay and its second moment (Ay)^. 

The calculation derived by Oide considers the beam 
size in the defocusing plane of the final quadrupole to 
be negligible, however, if the beam size in the defocusing 
plane is considered, the bending radius of curvature p in 
the lens is given by 

^ = |A:| ^x^s)+yHs) (9) 


and Fig. 3 shows schematically the particle trajectory. 
As shown in Eqs. (7) and (8), the radius of curvature 
depending on both transversal planes has a direct impact 
on the radiation. 

In addition, assuming I* /3*, the coordinates a;(s) 
and y(s) in Eq. (9) can be calculated from the horizon¬ 
tal {xq ) and vertical (j/q*) angles at the IP as 


cr(s)- 

(10) 

y{s)- 

(11) 


where 


fx{4>) = sinhcj) + Vlil* cosh (j) ( 12 ) 

Combining the expressions in the Eqs. (7), (9), (10) and 
(11), it is possible to average over the photon energy in 
Eq. (5), obtaining: 

A^ = a(y'*)3 + 5(a:'*)2y- (13) 

where 

2 py/kL -I 

a =3^67^ fyWFyi(l))d(j3= -rej^F^{VkL) (14) 

(y py/kL 

b=^rej^J fl{(j))Fy{4>)d(t) (15) 

and 

Fy^4>) = j%lW)dc^' 

-I- '/kl* sin^ (p. (16) 

Equation (13) shows a cubic component and in addition 
a linear component in j/q* whose magnitude depends on 
the horizontal angle at the IP, showing explicitly the cor¬ 
relation between the two planes when the beam size in 
the defocusing plane is not negligible. 

As we are interested in the effect over an ensemble of 
particles, the expected value of a function is defined as 

/ oo poo 

/ mx'o*,yo)dx'o*dy'o*, (17) 

-OO ^ —oo 

where 12 represents a Gaussian distribution in both hor¬ 
izontal and vertical particle angles at the IP. 

The average particle deviation due to radiation (Ay) 
is equal to zero, but in order to explore the correla¬ 
tion of the radiation with y^ the two components, a 
and 6(xo*)^, are evaluated further in this section using 
(( 2 ^ 0 *)^) = and later compared with 

particle tracking results of CLIC 3 TeV QDO. 

((Ay) ) is given by 

"“(I 

(18) 

and used in Section IIC to calculate the maximum theo¬ 
retical mitigation of the radiation effect and to evaluate 
the mitigation method. 

In a similar way, the square of the vertical dis¬ 
placement in Eq. (5), (Ay) , is calculated combining 
Eqs. (8), (9), (10) and (ll) by averaging over photons 
second moment of the energy loss, resulting in Eq. (19). 




















4 


( A ?/)2 


55 Ae 
24^3^® 27r 



(K/yW]'+ K/.('^)]') 


3/2 

Fim4>. 


(19) 


((Aj/)2) = when cr^;^. = 0. However, when is 
not zero the expected value cannot be calculated analyt¬ 
ically, it must be evaluated numerically and it is referred 
to as below. This is also compared with par¬ 

ticle tracking simulations in this section for the case of 
CLIC 3 TeV. 


B. Simulation 

Particle tracking from the entry of QDO to the IP 
for CLIC 3 TeV with and without radiation, using 
PLACET [10], allows to compute the effects of radia¬ 
tion on the six-dimentional phase space. Figure 4 shows 
the particle transverse distribution at the IP with and 
without radiation. 

40 
30 
20 

? 10 
0 

-10 
-20 
-30 

-0.2-0.15-0.1-0.05 0 0.05 0.1 0.15 0.2 0.25 

X (|tm) 

FIG. 4. CLIC 3 TeV beam at the IP after tracking through 
QDO with and without radiation. 

Figure 5 shows tracking results of Ay = y^ad — yno rad 
versus the particle angle yg at the IP for two different 
horizontal emittances. It has been included to ilustrate 
the radiation effect for the nominal normalized horizontal 
beam emittance, cpfx = 660 nm, and the case when the 
horizontal beam size Ux is negligible when compared to 
the vertical beam size which is achieved by reducing the 
horizontal emittance, tNx = 10“"^eAry, i.e. tNx = 2 pm. 
The horizontal beta twiss parameter at the IP is the same 
in both trackings, /I* = 6.9 mm. 

Fitting results of Ay for the two previously mentioned 
cases are compared in Table II with the analytical eval¬ 
uation of Eqs. (13), (14), (15) and (16). 

In the case ol cmx = 2 pm, i.e. a negligible effect 
of the defocusing plane, the result of a differs by 6% 
between theory and tracking results, while the result of 
biNx/ilPx) shows an unexpected component from track¬ 
ing. In the case oi e^x = 660 nm, i.e. when the defo- 
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FIG. 5. Correlation between the phase space coordinates 
Ay, y' for CLIC 3 TeV from particle tracking with two dif¬ 
ferent horizontal emittances. 
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a 

[10~“ m] 

beNx/i’fPx) 
[10"“ m] 

^Nx 

= 2 pm Theory 

9.0 

0 


Tracking 

9.5 ±0.1 

-1.3 ±0.3 

tNx 

= 660 nm Theory 

9.0 

6.3 


Tracking 

8.5 ±0.1 

5.4 ±0.3 


TABLE 11. Coeficients from fitting tracking results and theory 
Eq. (13): yo in [10“® rad] units. 


cusing plane is not negligible, the result of a differs by 
-6% between the theory and tracking results, while the 
result of b€Nx/{jl3x) shows -14% difference. In addition, 
there is an unexpected 11% variation in the fitting re¬ 
sults of the cubic component a caused by the change of 
the horizontal emittance in tracking. 

The contribution to the beam size ^(Ay)^^ from the 
two theoretical expressions (cToide and cr 2 D-oide), and the 
particle tracking are compared in Table III, where the er¬ 
rors included in the cr 2 D-oide results come from numerical 
precision in the integration of Eq. (19). 

The case of cnx = 2 pm shows that CToide and CT 2 D-oide 
agree within the numerical precision achieved. The re¬ 
sult from tracking is 6% above the theoretical values. 
On the other hand, the contribution to radiation calcu¬ 
lated from CT 2 D-oide with the nominal horizontal emittance 
£Nx = 660 nm agrees well with the tracking result. 

At the moment the differences between theory and 
tracking results for these cases have been attributed to 
limitations in the particle tracking and radiation Simula- 
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^oide 

0.87 


^ 2D-oide 

0.87 ±0.03 


Tracking 

0.92 

f-Nx = 660 nm 

^ 2D-oide 

1.02 ±0.03 
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1.00 


TABLE III. Contribution to beam size from radiation evalu¬ 
ated by (Joide, cr 2 D-oide and tracking. 



tions. 

Comparing (T 2 D-oide to CToide from Table III it is possi¬ 
ble to conclude that including the 2”*^ dimension has an 
important impact of 11% in the final IP vertical beam 
size in CLIC 3 TeV because of the additional 17% beam 
size contribution from radiation. 


C. Mitigating the impact on the beam size by a 
non-linear corrector scheme 

An ideal compensation system would remove the po¬ 
sition change due to radiation Ay = j/rad ~ yno rad- In 
this case, it is possible to remove the average effect due 
to energy loss Ay because it correlates with the particle 
angle in the vertical plane yg*. 

For CLIC 3 TeV, ((Ay)^)^^^ = 0.40 nm, composed 
of 0.11 nm from the linear and 0.34 nm from the cubic 
components in yg*. Using the expression 

V^jp-oide- ((^)") ^20) 

^2D-oide 

a possible ideal 8% reduction in the contribution to the 
Oide effect by removing the correlation. If only one of 
the two components is removed, then removing the linear 
component corresponds to 1% reduction, while removing 
the cubic components is 6% reduction in the contribution 
to beam size due to radiation. 

Several possibilities arise to achieve the mitigation of 
the beam size growth, in particular the beam size com¬ 
ponents due to fringe fields [11] show similarities with 
the components in Ay, pointing to the possibility to tune 
the lattice elements to compensate the effect of radiation, 
possibly using various octupoles. 

In this article, one possible correction scheme is tested 
consisting in the addition of one octupole magnet CO 
in front of the strong focusing quadrupole, placed as in 
Fig. 6. 

The kick given by the corrector, Ay'g, changes the 
vertical position at the IP as Ayco = Lc0^2/c0j where 
Lco is the distance from the IP to CO. The displacement 
induced by the corrector must have opposite direction to 


FIG. 6. The solid and dashed lines represent the particle tra¬ 
jectory without and with radiation, respectively. The dotted 
line represents the correction of the particle trajectory given 
by CO. 

the effect of radiation, Ayco = —Ay, therefore 

A2/'o = -^. (21) 

TcO 

The cubic and linear components of Ay suggest the use 
of an octupolar magnet as corrector, therefore, the kick 
given by an upright octupolar magnet is given by 


kz{x^ - ixy^) 

(22) 

ksi^x'^y-y^) 

(23) 


where fcg is the octupole strength. The difference in sign 
between the linear and cubic components when compar¬ 
ing Eq. (23) with Eq. (13) means that only one out of 
the two could be removed from the beam size by an oc¬ 
tupolar corrector, limiting the possible mitigation in this 
simple approach. 

The optimum location of the octupolar corrector is in 
front of the focusing magnet where the vertical beam 
size is large and the horizontal beam size is small which 
is achieved when maximizing the Py/jS^ ratio. Figure 7 
shows the horizontal and vertical /3 functions for the 
CLIC 3 TeV in the FD region, and their ratio, where 
= 5.3 X 10^ m and Py = 3.3 x 10® m at CO. 

Knowing the location of the corrector and matching 
the kicks in Eqs. (21) and (23) it is possible to obtain the 
two values for the octupole strength /cg, one to cancel the 
linear and another to cancel the cubic component. 

Linear: fcg = — (24) 

^cO 

f)n 

Cubic: fcg = (25) 

^cO 

where x'q = x/Lco and yg* = yjL^Q have been used. For 
CLIC 3 TeV and LcO = 3.5 m, the results of Eqs. (24) 
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FIG. 7. j3 functions and PyjPx ratio for CLIC 3 TeV FD, 
QFl and QDO in red on top. The dark area is occupied by 
QDO and the IP is at s = —3.5 m. 

and (25) is fca = — 2600 and = 3600 

respectively. 

In order to try this octupole mitigation in simulations, 
a CLIC 3 TeV nominal beam with no energy spread is 
generated at the IP and tracked back to the entrance of 
QDO without radiation with the corrector off. This beam 
is used to study the Oide effect mitigation in QDO using 
the previously mentioned scheme by tracking to the IP 
with radiation. The procedure consists in setting the best 
octupole strength for a 10 mm length CO. 

Results from the tracking to the IP with and without 
radiation with the corrector off, in Table IV, show that 
the Oide effect contribution to vertical beam size affects 
very little the total and peak luminosities, Ltot and Lp^aki 
obtained with Guinea Pig ++ [12]. 

The best result obtained with the octupole corrector is 
a vertical beam size reduction by (—4.3 ± 0.5)%, equiv¬ 
alent to 6% reduction in the Oide effect contribution to 
beam size which agrees with the removal of the cubic 
component on the beam size. 
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[nm] 

[nm] 
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NO RAD 

0 

47.45 

0.69 

7.7 

2.9 

RAD 

0 

47.45 

1.18 

7.5 

2.7 

RAD 

3900 

47.45 

1.13 

7.4 

2.7 


TABLE IV. Effect of the octupolar corrector on the beam 
size, total luminosity and peak luminosity. 

The CO strength from tracking is positive also indicat¬ 
ing the correction of the cubic component. However it is 
8% bigger than the theoretical value. This has been at¬ 
tributed to limitations in the simulation and the strength 
optimization. 

III. CONCLUSIONS 

Synchrotron radiation in the final quadrupole sets a 
lower limit in the beam size, which is particularly relevant 
for CLIC 3 TeV. This is normally mitigated by increas¬ 
ing the length of the magnet, nevertheless this produces 
larger chromatic and geometric aberrations. Previous re¬ 
sults considered only the focal plane of the quadrupole. 

A new value, cr 2 D-oide 5 has been derived to include 
the radiation effect of both the focusing and defocusing 
planes of the quad showing an increase of 17% in the con¬ 
tribution to beam size. Theory shows good agreement 
with simulations using the tracking code PLACET. 

A closer insight into the beam correlation induced by 
synchrotron radiation at the IP have revealed a way to 
mitigate the 2D-Oide effect with octupolar magnets. In 
general the 2D-Oide average geometric aberrations can 
be included in codes like MAPCLASS [6] to mitigate all 
aberrations simultaneously. The simple case of an oc¬ 
tupole is calculated and simulated giving a vertical beam 
size reduction of (4.3 ± 0.5)%, with little or negative im¬ 
pact on luminosity. 
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